Numerical study of fractionalization in an Easy-axis Kagome antiferromagnet 
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Based on exact numerical calculations, we show that the generalized Kagome spin model in the 
easy axis limit exhibits a spin liquid, topologically degenerate ground state over a broad range 
of phase space. We present an (to our knowledge the first) explicit calculation of the gap (and 
dispersion) of "vison" excitations, and exponentially decaying spin and vison 2-point correlators, 
hallmarks of deconfined, fractionalized and gapped spinous. The region of the spin liquid phase 
includes a point at which the model is equivalent to a Heisenberg model with purely two-spin 
interactions. Beyond this range, a negative "potential" term tunes a first order transition to a 
magnetic ordered state. The nature of the phase transition is also discussed in light of the low 
energy spectrum. These results greatly expand the results and range of a previous study of this 
model in the vicinity of an exactly soluble pointi 
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FIG. 1: Kagome lattice and interactions. Two primitive vec- 
tors a\,a2 are shown, as are the three 2-spin couplings on 
sites within a Hexagon. The ring term involves four sites on 
a bow-tie, which is generated from 2-spin virtual exchange 
processes. Sublattice labeling is also shown. 

The discovery of the fractional quantum Hall effect 
(FQHE) has revealed that electronic systems can have 
fractionally quantized excitations in nature. The possi- 
bility of such fractionalization in spin systems, without 
any applied magnetic field, is a subject of great inter- 
est. Anderson first proposed that two-dimensional (2d) 
spin-1/2 antiferromagnets might condense into a feature- 
less "spin-liquid" quantum ground stated, with decon- 
fined spinon excitations carrying spin S = 1/2^. In the 
past several years, a clear notion of fractionalized and 
gapped^ spin-liquid states emerged. With the absence 
of spin ordering and spatial symmetry breaking, such a 
liquid state is characterized by "topological order"—, as 
in the FQHEi^. Spinous are subject to long-range statis- 
tical interactions with vortex-like excitations (denoted as 
visons) which carry an Ising or Z2 fluJsSiiS. 

Theoretically a few spin models have been identified 
as possible candidates for realizing such a spin liquid 
pjjg^ggijiiji2ji3ji4^ Moessner and Sondhii^ have suggested 
this might occur for some antiferromagnets on the tri- 
angular lattice, by showing that a particular so-called 
"quantum dimer model" on the triangular lattice is in a 
featureless deconfined spin liquid phase in a range of pa- 



rameters. Unfortunately, the triangular quantum dimer 
model does not derive directly from a triangular lattice 
spin model, so which spin Hamiltonian might realize this 
state is currently unclear. Similar dimer models on the 
square lattice display only a deconfined critical point^ 
separating two confined phases. 

Balents, Fisher and Girvini have recently displayed 
a spin model on the Kagome lattice, closely mathe- 
matically related to the above triangular dimer model, 
which demonstrably has a fractionalized, topologically 
ordered ground state. Similarly to earlier studies on var- 
ious dimer modelsiSii^, fractionalization is established 
through mapping the system to an exactly soluble point, 
first exploited by Rokhsar and Kivclson (RK)4 in the 
square lattice. A drawback of this approach is that 
the RK point requires four-spin couplings, and topolog- 
ical order is argued perturbatively from the RK point. 
RefQ speculated, however, that the spin liquid state may 
persist to a simpler limit describable by only two-spin 
Heisenberg interactions. This central issue remains un- 
settled. 

In this letter, we present exact numerical diagonaliza- 
tion studies of a generalized Kagome spin modeli in the 
easy axis limit (see below) which interpolates between the 
RK point and a two-spin Heisenberg form - and beyond. 
We have found that the spin liquid phase persists over 
a wide range of parameters including the two-spin limit. 
This phase is characterized by the absence of spin order- 
ing, and by deconfined, fractionalized and gapped spinon 
excitations. The spin liquid phase has four-fold topolog- 
ical degeneracy, a finite gap to the vison excitations, and 
short range exponentially decaying spin and vison 2-point 
correlators. Still further from the RK point beyond the 
2-spin model, a first order transition to a magnetically 
ordered phase occurs. 

The model considered in most of this paper is the gen- 
eralized spin-1/2 ring-exchange Hamiltonian: 

H = ^{-JringS^ 828^ S;^ + h.C. + UiPfHp), (1) 
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where the labels 1 ... 4 denote the four spins at the ends 
of each bowtie (and others obtained by 120 degree ro- 
tations), as labeled in Fig. The additional term is a 
projection operator, 

PfHp^\nn){nn\ + \ini) ami, m 

with spin states iSf , 5"! , iSf , 5*1 indicated sequentially in 
the bras/kets. This Hamiltonian acts in the reduced 
Hilbert space with the constraint that for each hexagon, 
the total S*^ of six spins, Sq — 0. For the particular 
value U4 = 0, Eq. can be shown to be equivalent 
to the leading-order effective Hamiltonian describing the 
easy-axis limit of the Heisenberg model, 

n^2j2 Jt,si:s^, (3) 

where the sum is over pairs of sites (ij), with non-zero 
Jj^- = for all first, second, and third nearest-neighbor 
pairs (ij) on the Kagome lattice (see Fig. Specifi- 
cally, in the extreme easy-axis limit, Jz ^ J±, Eq. Q 
reduces by second order degenerate perturbation theory 
to Eq. ^ with Jring = ^Ji^l Jz ^-^^^ M4 = 0. The energy 
of states with 7^ is higher by 0( J^), and such states 
require additional terms beyond those in Eq. |^ for their 
description. 

The exact soluble RK point corresponds to M4 = Jring ■ 
For U4 = Jring — with e ^ Jring, the ground state is a 
featureless spin liquid state with gaps to all excitations. 
In particular, the vison gap was argued variationally to 
be 0{Jring)^- Wc havc performed exact Lanczos diago- 
nalization in the whole range of —2 < M4 < 1 (we take 
Jring = 1 as the Unit). Specifically, we consider a finite 
size system on the torus with length vectors Li = nicii 
and L2 = ?^2a2l which connect identical sites (periodic 
boundary condition (PBC)). Here di and 0,2 are the prim- 
itive vectors shown in Fig. 1, and we set lattice constant 
fli = 02 = 2 for convenience. Then total number of sites 
is Ns = 3niri2. In the constrained Hilbert space with 
Sq = 0, this problem is characterized by two topological 
Z2 "winding numbers" , (101,^2) 

Wa^^2Sl=±l, (4) 
k 

where, for concreteness, the product is taken on one ar- 
bitrarily chosen straight line of sites along the ai /a2 axis 
for a — 1,2 encircling the torus. These two winding 
numbers are complete in that other winding numbers de- 
fined on other non-trivial loops are not independent of 
these. Using translational invariance, the Hilbert space 
for Ns = 60 can be reduced to a dimension (total number 
of configurations) around Nconf — 3618896 {Nconf varies 
amongst different topological sectors). 

Shown in Fig. 2 is the excitation energy gap between 
the lowest two states Eg — E{2) — E{1) in the topological 
sector with winding numbers (1,1)- As U4 moves away 
from 1, the energy gap remains at the order of 1 and 
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FIG. 2: The spectral gap Eg = E{2) - E{1) in the winding 
sector (1,1) as a function of U4 for sizes Ns = 36,42,48,54 
and 60. In the inset, AEg/Au4, vs. U4 for Ns =42 — 60 (same 
symbols as in the main plot). 

only becomes much smaller on the negative U4 side. All 
curves actually cross over each other around Uc — —0.5, 
and Eg overall decreases with the increasing of Ng with a 
trend of going to zero in the regime U4 < Uc = —0.5. To 
examine the finite-size effect of transition, we show the 
slope of the Eg curve vs. U4 in the inset of Fig. 2 for large 
sizes Ns = 42 — 60. Indeed around Uc ~ —0.5, a strong 
peak showed up in AEg/ Au4, with its height increasing 
with Ns, which is suggestive towards a quantum phase 
transition at Uc to a different phase with vanishing Eg. 

At the RK point, a finite gap Eg is expected, related 
to the excitation energy of vison excitationsi, and as just 
discussed, this gap persists for Uc < U4 < Urk = 1- 
Visons are characterized by the Z2 flux $^ = ±1, which 
is defined on each triangular plaquette of the Kagome 
lattice, asi 

'^''^ = n n 2^^^ (5) 

where the first product is over the three bowties centered 
on the sites of the triangle, and the second is over the 
sites of each bowtie. It is simple to show that a ground 
state |0) of H can always be found which is expressed 
as a superposition of S^ eigcnstatcs with positive real 
coefficients. This implies that (0|$a|0) > 0; hence, there 
are no visons in the ground state. One can readily show 
that with PBCs the product of <I>a over all triangular 
plaquettes (even over just all say up-pointing triangles) 
is -f-l, hence visons (^a = — 1) can only appear in pairs 
with PBCs for Ti,. An appropriate definition of a single 
vison state is made as follows. We imagine a large open 
system, and perform the canonical transformation 

H' = VioHvio, (6) 

where Vio denotes a single vison creation operator iden- 
tified in Ref. Q, which is a "string" operator made of the 
product of spins 2S'f along some path on the kagome 
lattice starting at site iq and ending at the boundary. 
Explicitly, H' is the same as the TC except that the ring 
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FIG. 3: (a) The single vison energy Esv = £'(0) - £'(0) as a 
function of Ui for system sizes Ns = 36, 42, 48, 54 and 60; (b) 
The low energy spectrum Ai5„ = E{n) — E ground for Ns = 60 
(Li = 5ai and L2 — 4a2), in all four sectors vs. U4. 



exchange terms of the three bow-ties centered on the tri- 
angle containing iq from which the path exits are changed 
in sign. This canonical transformation redefines the Z2 
flux on this triangle only by a minus sign. $a — > — ^A, 
i.e. a vison on this triangle corresponds now to $a = +1. 
Imposing PBCs on TL' (it is no longer canonically con- 
jugate to Ti) then forces the total number of visons to 
be odd. Although it is not manifest, TC' has a hidden 
"magnetic" translational invariance, so there is no lo- 
calization of the forced vison. Thus the energy differ- 
ence between the ground states of the two Hamiltonians, 
Esv = E'{0) — E{0) gives the single vison energy. We 
plot Esv vs. U4 in Fig. 3a for system sizes A^g = 36 — 60. 

Esv shows overall similar behavior to the spectral gap 
Eg, but bigger than 0.5Eg in the vison gapped regime, 
indicating a non-vanishing binding energy between two 
visons. In the vison gapped regime, we find that the low 
energy manifold of Ti.' forms an energy band of single vi- 
sons with very small dispersion (eg., for Ns = 60 disper- 
sion of single vison energy is about 0.047 Jring at U4 = 0). 
Esv remains to almost constant in the regime passing 
M4 = 0, and drops significantly at 114 < J7c = — 0.5 sidsi^. 

To reveal the nature of the phase transition at Uc, we 
present the low energy spectrum for Ns = 60 in all four 
topological sectors by plotting AEn = E{n) — E ground 
vs. U4 in Fig. 3b (where Eground refers to the lowest 
state energy of the system). Clearly, in the whole range 
of U4 > Uc, the lowest energy states from each topolog- 
ical sector remain quasidegenerate (filled triangular-up, 
is the energy difference). At U4 < Uc, the states with dif- 
ferent winding numbers along Li, are separated into two 
groups. All the low energy states from the sectors (1, ±1) 
go down and collapse with the ground state. The other 
(—1, ±1) states rise and become well separated from the 
low energy manifold. This dependence on the topological 
sector is evidence of developing spin long-range ordering. 
Similar results are obtained for all system sizes 36 — 66. 
The phase transition between the vison gapped state and 
spin ordered phase seems likely first order as Eg crosses 
at Uc and becomes zero at U4 < Uc side. At almost 
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The spin correlator dj vs. distance Xi 
—0.4, and —1 for strip- like system at Ns — 
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the same time, the ground state degeneracy of all sectors 
breaks done and long-range spin order develops. 

As discussed in Ref. |^ exponential decay of the vison- 
vison correlation function is the hallmark of a 2d Z2 frac- 
tionalized phased. Following Ref. we define the spin 
and the vison 2-point correlation functions: 



lH>25£:|0) 



(7) 



where |0) denotes the ground state, and the product in 
Vij is taken along some path on the Kagome lattice start- 
ing at site i and ending at site j, containing an even 
number of sites, and making only "±60°" turns left or 
right. Due to the constraint = 0, the latter prod- 
uct is path-independent up to an overall sign (hence the 
absolute value in Eq. [Tj) . 

To examine the longer distance behavior, we consid- 
ered a strip geometry with L2 = 2a2, while varying 
Li ~ nidi. Similar analysis of the spectral gap Eg and 
single vison energy Esv reveals that the critical Uc for 
strip-like samples is slightly more negative than that for 
more two-dimensional clusters, around —0.75. 

The numerically calculated Cij is shown in Fig. 4. for 
Ns = 66 and different U4. At U4 = 1 (RK-point), Cij 
clearly shows the exponential decay seen previously from 
the exact wavefunctioni. As U4 varies from 1 to 0, apart 
from a small oscillation, we see essentially the same ex- 
ponential behavior. The data at M4 = 1 and 0, can both 
be well fitted by InCi^ ~ — jxi — with apparently 

the same correlation length ^ w 1.7. For U4 further de- 
creased to —0.4, just before the phase transition, much 
stronger fluctuations emerge between Xi — Xj even or odd, 
but the overall decay remains exponential. At U4 = —1, 
Cij jumps by more than one order of magnitude at longer 
distance Xi — Xj ~ Li/2 = 11, and a longer range corre- 
lation is clearly evident. 

We analyze the finite-size dependence of the spin spin 
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FIG. 5: (a) dj at fixed Xi — Xj = L\/2 vs. it4 for strip- 
Uke system at Ns ~ 36, 42, 48, 60, 66. (b) Vison correlator Vij 
along ai (vison-1) and 02 (vison-2) for 2d- like system Ns = 60 
are compared with dj for both 2d system at Ns = 60 and 
stripe system at A'^s ~ 66. 

correlation function. In Fig. 5a, Cy with j = i + Li/2 
(halfway across the torus) is shown for relatively large 
sizes Ns = 36, 42, 48, 54, 60 and 66. For U4 > Uc (around 
-0.75), Cij is vanishingly small for all sizes. For U4 < Uc, 
however, Cij is very weakly dependent on Ns (except 
for the smallest system with Ns = 36), and it robustly 
scales to a finite value at large Ng limit. This strongly 
indicates long-range magnetic order. One expects this 
region is adiabatically connected to the U4 — > — 00 limit, 
for which the ground state is fully magnetically ordered 
state with 4:{SiSj) = ±1, taking the positive or negative 
sign if i and j belong to the same or different sublattices, 
respectively (see sublattice A and B labeling in Fig. 1; 
spins not belonging to these two sublattices are not or- 
dered). In the strip geometry for large negative 144, one 
can show the ground state manifold (with 4-fold degen- 
eracy) is well separated from excited states by an energy 
gap of 4Jrmg/|u4p, with additional ordering in the x- 
y plane (in spin space), as a result of order- by-disorder 
phenomena. Details will be reported elsewhere. 



We further examine the vison 2-point correlator. In 
the strip case, the vison correlator vanishes exponentially 
but with an extremely short correlation length ^ much 
smaller than the lattice constant. For more 2d clusters, 
the vison correlator has almost the same f as the spin 
correlator. We plot Cij and Vij at 1*4 = for cluster 
of Ns = 60, with Li = 5ai and L2 = 4a2, vs. Ri — 
Rj (along both ai or 0,2 as vison-1 or vison-2) in Fig. 
5b. Despite the very short distance across the torus, 
the several points indeed follow an exponential behavior 
comparable with Cij for the strip with Ng = 66. Clearly, 
all the correlators can be well fitted by exp(— ji?; — i?j|/^) 
(as the dashed line in the plot) with ^ k, 1.7. 

We conclude with some discussion of related quantum 
phase transitions. Our numerical calculations indicate 
that the fractionalized spin-liquid state, previously es- 
tablished for the generalized Kagome spin model at the 
RK point, persists in a wide range of parameters, cov- 
ering the pure 2-spin Kagome model (1/4 = 0). On the 
negative M4 side, a first order transition to the magnetic 
ordered phase happens at Uc = — 0.5 J^mg- It seems that 
the easy axis spin model easily develops long-range Ising 
magnetic order at the confinement transition. It is possi- 
ble that this ordering is an artifact of the easy- axis limit, 
which suppresses some fluctuations, and perhaps Kagome 
antiferromagnets away from easy-axis limit may develop 
instead a gapless spin-liquid phase, which would be ex- 
tremely interesting to explore. 
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